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Abstract 

We point out that applying the photon energy cut significantly modifies the moments of 
energy spectrum in B—>-Xs + '-f decays, with a certain class of effects not accounted for in 
the mostly used OPE expressions. This leads to a systematic bias in the extracted values of 
the b quark mass and other heavy quark parameters. The apparent b quark mass increases 
typically by 70 MeV or more, together with an even more dramatic downward shift in 
the kinetic expectation value. Accounting for these cut-related shifts brings different 
measurements into a good agreement, when the OPE-based theory employs the robust 
approach. These nonperturbative effects are exponential in the effective hardness severely 
lowered by high cuts, and do not signify a breakdown of the 1 / expansion itself. Similar 
effects in semileptonic — > c decays are briefly addressed. We stress the utility of the second 
moment of once these effects are incorporated. 
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Treating effects of strong interactions in decays of beauty particles is of primary impor- 
tance nowadays when many precision experimental studies are being carried out. Inclusive 
distributions in radiative and semileptonic decays are the portal to accurately determin- 
ing the nonperturbative heavy quark parameters controlling short-distance observables in 
B decays. Their utility rests on an consistent expansion in l/rrib, the inverse b quark 
mass [1] . Recent experimental data generally show a nontrivial agreement between quite 
different - and a priori unrelated - measurements at the nonperturbative level, on one 
hand, and consistency with the QCD-based OPE theory. 

In order to enjoy in full the potential of a small expansion parameter provided by 
the heavy quark mass, the observable in question must be sufficiently inclusive. However 
experimental cuts imposed for practical reasons - to suppress backgrounds etc. - often 
essentially degrade the effective hardness Q of the process. This brings in another expan- 
sion parameter 1/Q effectively replacing l/m^ in certain QCD effects. The reliability of 
the expansion greatly deteriorates for Q<^mb. This phenomenon is particularly impor- 
tant in 6 — > s -|- 7 decays where experiments so far have imposed cut E^>2 GeV or even 
higher. 

The theoretical aspects of such limitations have been discussed during the last couple 
of years [2, 3, 4]. In particular, the effective 'hardness' of the inclusive 6 — > s -|- 7 decays 
with Ej>2 GeV amounts to only about 1.25 GeV, which casts doubts on the precision of 
the routinely used expressions incorporated into the fits of heavy quark parameters. 

In the present letter we point out that these effects are numerically significant, lead 
to a systematic bias that often exceeds naive error estimates and therefore cannot be 
ignored. Evaluating them in the most straightforward (although somewhat simplified) 
way we find, for instance for 6 — > s -|- 7 decays 

rfib ~ mb + 70 MeV 

ill ^ /.^ - (0.15 - 0.2) GeV^ (1) 

where rhb and /i^ are the apparent values of the b quark mass and of the kinetic expec- 
tation value, respectively, as extracted from the 6 — > s -|- 7 spectrum with > 2 GeV 
in a usual way. Correcting for these effects eliminates alleged problems for the OPE in 
describing different data and rather leads to a too good agreement between the data on 
different types of inclusive decays. 

Moreover, this resolves the controversy noted previously: while the values of A and //^ 
reportedly extracted from the CLEO 6 — > s -|- 7 spectrum were found to be significantly 
below the theoretical expectations, the theoretically obtained spectrum itself turned out 
to yield a good description of the observed spectrum when we evaluated it based on these 
theoretically preferred values of parameters [2] . 

The bias in Eqs. (1) depends on the position of the cut (more precisely, on the gap 
£^cut) s-nd the actual values of other heavy quark parameters. The quoted estimates 
assume moderate values, mb{l GeV) ~4.6 GeV and /U^(l GeV) ~0.43 GeV^. If nib becomes 
lower and/or the true /^^ increases as may be indicated by the most recent data, the bias 
increases further. 
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OPE and cuts The origin of these effects and why they are missed in the standard 
apphcation of the OPE and in estimates of the theoretical accuracy, have been discussed 
elsewhere [4]. In brief, considering a constrained fraction of the B ^ Xg+'y events 



(or similarly truncated photon energy moments), the simple-minded approach routinely 
expands the spectrum in powers of l/nii. Ignoring perturbative bremsstrahlung one 
obtains a (5-like spectrum peculiar for two-body decays, and the expansion around the free- 
quark kinematics does not change this - it only generates higher derivatives of 5(£'-y— ^) : 

Tt,^ = «^(^7-f^) + b6'{E,-^) + c6"{E,-^-f) + ... (3) 

where a, b, ... are given by the B meson expectation values of local 6-quark operators. 
Naively computing 1— or spectral moments over the restricted domain in this way 
would yield unity in Eq. (2) for any E > ^ - a result clearly meaningless on physical 
grounds. The actual behavior of the spectrum and the moments is described by the heavy 
quark distribution function. Its tail is indeed exponentially suppressed by a typical factor 
g-cQ/Mhadr Q(-B-y) ^ A*hadr ■ However, for ^(-B^) ~ ^hadr there is little suppression of 
the missed tail contribution; the error becomes of order one. 

This obvious point is missed in the naive application of the OPE and in the way to 
gauge the theoretical uncertainty based on it. Conceptually this is related to the limited 
range of convergence of the OPE for the width, determined in this case by the support of 
the heavy quark distribution function [4]. 

Here we rather concentrate on the numerical consequences for B^Xs + ^ decays. To 
this end we first turn off perturbative effects altogether. The spectrum then is described 
by the nonperturbative light-cone distribution function F{k^): 

^^ = 2F(2i5,-m,). (4) 

Although not necessary for our purpose, one can imagine a theoretical heavy quark limit 
with fixed hardness Q : 

Q = nib — 2Ecut = fixed , mf, — >• oo . (5) 

The fully integrated moments of F{x) and therefore of the spectrum give then directly 
the underlying heavy quark parameters: 



/•oo roo /|2 

/ kF(A-k)dk = A, / (k-A?F(A-k)dk = ^, etc. (6) 

^0 Jo 3 

As mentioned above, in the standard practical-OPE-based formulae these relations remain 
the same for the moments with the cut as well (provided E^ut < ^ which is always 
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assumed) - yet not in reality. Paralleling the routinely used way we therefore introduce 



cut; 



\2 dr 



dE. 



-A'(i5^cut) 



(7) 

Clearly the apparent value A(i?cut) is always below the actual A. This is illustrated by 
Fig. la where we plot the cut-related 'bias' - the difference between A and K{E 

cut / 3j 

function of energy E^ut- It turns out to be quite significant. 
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introduced by imposing a lower cut in the photon energy in i? ^ Xg+'j. Blue and maroon 
curves correspond to two different ansatze for the heavy quark distribution function, Fi and F2, 
respectively. 



The naive extraction of the kinetic expectation value through the variance of the 
truncated distribution undercounts it even more dramatically as Fig. lb illustrates, since 
higher moments are more sensitive to the tail of the distribution. 

To gauge the sensitivity to the choice of the heavy quark distribution function we 
follow Ref. [2] evaluating the effect for two ansatze - one exponential in k^, and the other 
in k"^ yielding an even faster decreasing tail: 

Fi(A;+) = iVi (A-k+re''+ 9{K-k+), F2{k+) = N2 (K-k+f e-'^'^'^-^^'^" e(A-k+) ; (8) 

the parameters are adjusted in such a way as to yield the same actual rrih (A) and /i^. 
(For mfe = 4.6 GeV and /i^ = 0.43 GeV^ we have a = 2.22 and /3 = 0.773 .) Curiously, the 
'deficit' A— A(£'cut) practically does not depend on the choice at Ecut around 2 GeV, and 
even the deficit in /i^ is reasonably stable. 

Why are the effects of the cut so significant? They are exponential in the inverse 
hadronic scale /Xhadr, but are governed by the hardness Q~mb— 2ii^cut rather than by m;,: 

_ Q Q 

A-A(Ecut) OC /ihadr e"'^hadr , nl-fj,l (X /Uhadr e'^hadr (9) 
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(the exponent may have a form of a power of Q//Xhadr)- Even at mi, — > oo these effects 
survive unless Q is made large as well! 

As explained in Ref. [4] the bias terms (9) are associated from a theoretical viewpoint 
with the limited (in fact, zero) convergence radius of the OPE. This becomes practically 
relevant due to the presence of a subseries in powers of 1/Q rather than l/mi,. The 
limitations on convergence appear due to a factorial growth of the power coefficients, a 
rather universal property of the OPE. In this respect one may associate this effect with 
quark-hadron duality [5]. Yet it has no features peculiar to local quark hadron duality 
violation intrinsic to inclusive decay widths in actual Minkowski world. (For a discussion 
of the notorious subtleties in the notion of quark-hadron duality, see reviews [6, 5]). For 
instance, these effects are truly exponential and do not oscillate.^ 

The validity of the routinely applied expressions for the moments with cuts for actual 
decays is additionally complicated by perturbative corrections. Incorporated into fits are 
naive sums of pure perturbative and pure nonperturbative terms: 



where nonperturbative corrections to the moments M„ still do not depend on Ecut- This 
is not true in general, but would hold if the actual spectrum were exactly a convolution 
of the perturbative and nonperturbative spectra. 



provided no cut is introduced (or its effect on the pure nonperturbative distribution is 
negligible). We hasten to add, though that the effects we consider are unrelated to this 
complication and rather represent an independent phenomenon - they are significant even 
in the complete absence of perturbative corrections. 

Since the perturbative effects can potentially modify the effect, we have evaluated 
the cut-induced deficit in A and fl^ including short-distance corrections. Namely, we 
considered A, /i^ in Eqs. (7) for the complete spectrum obtained by the convolution (11) 
of the perturbative and primordial (nonperturbative) ones, and compared them to the 
naive sum Eq. (10) which would indeed hold for a sufficiently low cut. Including the 
perturbative spectrum as detailed in Ref. [2] we found no appreciable change in A— A or 
//^ — /t^ at realistic cuts (a small increase emerged only at £^cut;^l GeV). 

Based on these results we conchide: 

• the value of rrib as routinely extracted from the b^s + j spectrum is to be decreased 
by an amount of order 70 MeV; 

• relative corrections to are even more significant and can naturally constitute a shift 
of 0.2 GeV^. This arises on top of other potential effects. 



Practical implications Accepting the above shifts at face value and using the rather 
arbitrary choice = 4.595 GeV, = 1.15 GeV, //^ = 0.45 GeV^ = 0.06 GeV^ and 



^Peculiarities of real local duality would appear here only at the next-to- leading order in 1/mb and, 
therefore are not of much interest. 




(10) 




(11) 
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p\g = —0.15GeV'^ adjusted to accommodate {M]^) Eg^iGeV , we obtain^ 

(M^) ^ 4.434 GeV^ [cf. (4.538 ± 0.093) GeV (New DELPHI) ] 

{M]^)E,>i.5GeV ^ 4.177 GeV^ [cf. 4.180 GeV^ (BaBar), 4.189 GeV^ (CLEO) ] 

{Ee) ~ 1.389 GeV [cf. (1.383 ± 0.015) GeV (DELPHI) ] 

{E^)E,>2GeY ^ 2.329 GeV [cf. (2.346 ± 0.034) GeV (CLEO) ] 

{E^-E^)E,>2GeY ^ omll GeV^ [cf. (0.0226 ± 0.0066 ± 0.0020) GeV' (CLEO) ] (12) 

(experimental data are from Refs. [7, 8, 9, 10]). ii^f^t-dependence of {Mx) is also repro- 
duced, see Fig. 2. The counterpart of the above exponential cut-related effects for the 
semileptonic decays has not been incorporated here, however. 




Figure 2: Experimental values of the average hadronic mass square (M|.) in B ^ Xciu at 
different lower cuts in and the literal OPE prediction (blue curve) for the stated heavy quark 
parameters. The DELPHI point assumes no cut on E^. 

It should be noted that the cut-induced shifts rhb—rrib, /i^— /^^ are not unambiguously 
determined in terms of a few known HQ parameters, but rather depend on the actual shape 
of the heavy quark distribution function; in particular it is driven by the asymptotics of 
certain high- dimension expectation values. The evaluation presented above relies on the 

^The two values for the second iS^-moment correspond to the two ansatze; they are obtained discarding 
higher-order power corrections to the light-cone distribution function. 
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most natural assumptions about the function. Strictly speaking, they can vary for more 
contrived ansatze. Therefore, the estimates in Eqs. (1) can be conservatively viewed as 
the minimal possible theoretical inaccuracy of the usual naive evaluations. The safest way 
to tackle possible ambiguity is to resort to more inclusive moments without too severe 
cuts. 

Cuts in semileptonic moments. As argued in Ref. [4] similar cut-related 'exponen- 
tial' biases missed in the naive OPE applications affect the truncated moments in the 
semileptonic decays as well. Their description, even simplified is less transparent and 
would be more involved, though. In particular, the light-cone distribution functions for 
6— >s -|- 7 is replaced by a different function, the form of which actually is not universal. 
Yet the qualitative trend is expected the same - it should smoothly interpolate the case 
of literal OPE expression at low cuts and the values at a high cut dictated simply by 
the actual hadrons kinematics. This would replace step- or (5-like behavior in the formal 
OPE expressions. 

Numerical aspects are less certain. Keeping in mind that for B — > the effective 
hardness 



at E'cut = 1-5 GeV is about 1.25 GeV [3] , nearly the same as for S — > X5-I-7 with ^2 GeV, 
we may expect quite significant effects. To get a rough idea of the possible magnitude 
of the bias we can use a simplified rule of thumb - assume that nib in the semileptonic 
decay can be just replaced by an effective larger value rfife apparent in 6 — > s -|- 7 at the 
commensurate cut yielding the same hardness. In other words, we mimic the effect of 
decreasing hardness by an additional effective nonperturbative running of the heavy quark 
mass at low scales. As illustrated above, this is to increase m;, (or, equivalently decrease 
A) by about 70 MeV, a significant change. 

Yet the semileptonic decay characteristics strongly depend on both and itlc. To 
stay on the conservative side we assume the apparent shift in rric as high as in (which 
would reflect just the heavy flavor symmetry): 



In actuality the corrections to rric are typically somewhat smaller due to 1/m^ terms. 

As discussed elsewhere [11], for actual B decays both lepton moments and {Mx) 
depend on more or less the same combination of masses m^ — CGSm^. This means that 
the literal ansatz (14) would suppress the effect by a factor of 3 to 4, yet l/rric corrections 
may be thought to softening this suppression. We then expect the exponential terms in 
semileptonic decays with E'cut — l-5GeV to introduce effects on the same scale as shifting 
rrib upward by up to 25 to 30 MeV (assuming fixed rric and other heavy quark parameters). 
This rule of thumb is useful to get an idea of the ultimate theoretical accuracy one can 
count on. 



For example, the CLEO's cut moment Ri = {Ee)^^^^^^ is approximately given by [3] 
= 1.776 GeV + 0.27(m6-4.595 GeV) -0.17(mc-l. 15 GeV) at iKb/K&HO.OS (15) 




(13) 



rric TTT'c ^ mc+ {mb—mb) ■ 



(14) 
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(the above mentioned values of the nonperturbative parameters are assumed). An increase 
in ruf) by only 20 MeV would then change 

i?i ^ i?i + 0.0055 GeV (16) 

and would perfectly fit the central CLEO's value 1.7810 GeV. It is worth noting that 
Eqs. (15)-(16) make it explicit that the imposed cut on degrades the theoretical calcu- 
lability of Ri far beyond its experimental error bars, the fact repeatedly emphasized over 
the last year. Unfortunately, this was not reflected in the fits of parameters which placed 
much weight on the values of Rq - R2 just owing to their small experimental uncertainties, 
whilst paying less attention to actual theoretical errors.^ 

Similar reservations apply to the quality of theory for CLEO's R2 representing the 
second moment with the cut, with the only difference that the effective hardness only 
deteriorates for higher moments. 

The ratio Rq is the normalized decay rate with the cut on Ei as high as 1.7 GeV, and 
for it hardness Q is below 1 GcV. A precision - beyond just semiquantitative - treatment 
of nonperturbative effects is then questionable, and far more significant corrections should 
be allowed for. 

We note that there is a good agreement of most data referring to sufficiently 'hard' 
decay distributions with the theory based on the OPE in QCD, in the 'robust' OPE 
approach. The latter was put forward [3] to get rid of vulnerable and unnecessary as- 
sumptions of the usually employed fits of the data. This consistency likewise refers to 
the absolute values of the heavy quark parameters necessary to accommodate the data. 
Theory anticipates, however that the expansion becomes deceptive with increase of the 
experimental cuts. Here we have addressed the most obvious effects, those from the va- 
riety of 'exponential' terms in the effective hardness. While presently not amenable to 
precise theoretical treatment, they can be estimated using quite natural assumptions and 
are found to be far too significant for E^^^ ^1.5 GeV and E^^t ^2 GeV often employed 
in experiment. Taking these corrections at face value and incorporating in our predic- 
tions, we get a good, more than qualitative agreement with "less short-distance" inclusive 
decays as well. 

To summarize, the cuts essentially decreasing the hardness in B decays introduce 
terms some of which are exponentially suppressed though only in the effective hardness, 
but not suppressed by powers of the heavy quark mass. They significantly change the 
extracted values of the heavy quark parameters leading to an apparent suppression of 
the magnitude of the nonperturbative ones in B mesons. Accounting for such effects 
appears necessary in B ^ Xg + ^ decays unless the cut on photon energy is pushed 
well below present 2 GeV. This brings the practical predictions for various inclusive 
characteristics into a good agreement. Moreover, it seems that the second moment of 
the photon spectrum discarded for the fits so far, does impose informative constraints 

^N.U. acknowledges that a similar in spirit criticism of the theory error treatment in [12] was expressed 
by D. Hitlin at the BaBar Workshop, SLAG December 2002. 



7 



provided these effects are properly incorporated. 
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